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Abstract 

We study the problem of subharmonic bifurcations for analytic systems in the plane with pertur- 
bations depending periodically on time, in the case in which we only assume that the subharmonic 
Melnikov function has at least one zero. If the order of zero is odd, then there is always at least one 
subharmonic solution, whereas if the order is even in general other conditions have to be assumed to 
guarantee the existence of subharmonic solutions. Even when such solutions exist, in general they 
are not analytic in the perturbation parameter. We show that they are analytic in a fractional power 
of the perturbation parameter. To obtain a fully constructive algorithm which allows us not only to 
prove existence but also to obtain bounds on the radius of analyticity and to approximate the solu- 
tions within any fixed accuracy, we need further assumptions. The method we use to construct the 
solution - when this is possible - is based on a combination of the Newton-Puiseux algorithm and the 
tree formalism. This leads to a graphical representation of the solution in terms of diagrams. Finally, 
if the subharmonic Melnikov function is identically zero, we show that it is possible to introduce 
higher order generalisations, for which the same kind of analysis can be carried out. 

1 Introduction 

The problem of subharmonic bifurcations was first considered by Melnikov [S], who showed that the 
existence of subharmonic solutions is related to the zeroes of a suitable function, nowadays called the 
subharmonic Melnikov function. The standard Melnikov theory usually studies the case in which the 
subharmonic Melnikov function has a simple (i.e. first order) zero In such a case the problem can 

be reduced to a problem of implicit function theorem. 

Nonetheless, it can happen that the subharmonic Melnikov function either vanishes identically or has 
a zero which is of order higher than one. In the first case hopefully one can go to higher orders, and if 
a suitable higher order generalisation of the subharmonic Melnikov function has a first order zero, then 
one can proceed very closely to the standard case, and existence of analytic subharmonic solutions is 
obtained. Most of the papers in the literature consider this kind of generalisations of Melnikov's theory, 
and often a second order analysis is enough to settle the problem. 

The second case is more subtle. The problem can be still reduced to an implicit function problem, 
but the fact that the zeroes are no longer simple prevents us from applying the implicit function theorem. 
Thus, other arguments must be used, based on the Weierstrass preparation theorem and on the theory of 
the Puiseux series [HI El |3l H] . However, a systematic analysis is missing in the literature. Furthermore, 
in general, these arguments are not constructive: if on the one hand they allow to prove (in certain cases) 
the existence of at least one subharmonic solution, on the other hand the problem of how many such 
solutions really exist and how they can be explicitly constructed has not been discussed in full generality. 



The main difficulty for a constructive approach is that the solution of the implicit function equation 
has to be looked for by successive approximations. At each iteration step, in order to find the correction 
to the approximate solution found at the previous one, one has to solve a new implicit function equation, 
which, in principle, still admits multiple roots. So, as far as the roots of the equations are not simple, 
one cannot give an algorithm to produce systematically the corrections at the subsequent steps. 

A careful discussion of a problem of the same kind can be found in [T] , where the problem of bifurca- 
tions from multiple limit cycles is considered - cf. also [Ulin], where the problem is further investigated. 
There, under the hypothesis that a simple (real) zero is obtained at the first iteration step, it is proved 
that the bifurcating solutions can be expanded as fractional series (Puiseux series) of the perturbation 
parameter. The method to compute the coefficients of the series is based on the use of Newton's polygon 
[H [31 [T], and allows one to go to arbitrarily high orders. However, the convergence of the series, and 
hence of the algorithm, relies on abstract arguments of algebraic and geometric theory. 

To the best of our knowledge, the case of subharmonic bifurcations was not discussed in the literature. 
Of course, in principle one can think to adapt the same strategy as in T for the bifurcations of limit 
cycles. But still, there are issues which have not been discussed there. Moreover we have a twofold aim. 
We are interested in results which are both general - not generic - and constructive. This means that we 
are interested in problems such as the following one: which are the weaker conditions to impose on the 
perturbation, for a given integrable system and a given periodic solution, in order to prove the existence 
of subharmonic solutions? Of course the ideal result would be to have no restriction at all. At the same 
time, we are also interested in explicitly construct such solutions, within any prefixed accuracy. 

The problem of subharmonic solutions in the case of multiple zeroes of the Melnikov functions has 
been considered in [H], where the following theorem is stated (without giving the proof) for C" smooth 
systems: if the subharmonic Melnikov function has a zero of order n < r, then there is at least one 
subharmonic solution. In any case the analyticity properties of the solutions are not discussed. In 
particular the subharmonic solution is found as a function of two parameters - the perturbation parameter 
and the initial phase of the solution to be continued -, but the relation between the two parameters is 
not discussed. We note that, in the analytic setting, it is exactly this relation which produces the lack 
of analyticity in the perturbation parameter. Furthermore, in [14 the case of zeroes of even order is not 
considered: as we shall see, in that case the existence of subharmonic solutions can not be proved in 
general, but it can be obtained under extra assumptions. 

In the remaining part of this section, we give a more detailed account of our results. One can formulate 
the problem both in the C" Whitney topology and in the real-analytic setting. We shall choose the latter. 
From a technical point of view, this is mandatory since our techniques requires for the systems to be 
analytic. However, it is also very natural from a physical point of view, because in practice in any physical 
applications the functions appearing in the equations are analytic (often even polynomials), and when 
they are not analytic they are not even smooth. Also, we note since now that, even though we restrict 
our analysis to the analytic setting, this does not mean at all that we can not deal with problems where 
non-analytic phenomena arise. The very case discussed in this paper provides a counterexample. 

We shall consider systems which can be viewed as perturbations of integrable systems, with the 
perturbation which depends periodically in time. We shall use coordinates (a. A) such that, in the 
absence of the perturbation, A is fixed to a constant value, while a rotates on the circle: hence all 
motions are periodic. As usual [7] we assume that, for A varying in a finite interval, the periods change 
monotonically. Then we can write the equations of motion as d = iij{A) + eF{a, A, t), A — eG{a, A, t), 
with G, F periodic in a and t. All functions are assumed to be analytic. More formal definitions will be 
given in Section [2l 

Given a unperturbed periodic orbit t —^ (ao(i), ^o(i))i we define the subharmonic Melnikov function 
M{to) as the average over a period of the function G{ao{t), Ao,t + to). By construction Mito) is periodic 
in to- With the terminology introduced above, e is the perturbation parameter and is the initial phase. 



2 



The following scenario arises. 

• If Af(to) has no zero, then there is no subharmonic solution, that is no periodic solution which 
continues the unperturbed one at e 7^ 0. 

• Otherwise, if M{to) has zeroes, the following two cases are possible: either M{to) has a zero of 
finite order n or M(to) vanishes with all its derivatives. In the second case, because of analyticity, 
the function M{to) is identically zero. 

• If M{to) has a simple zero (i.e. n = 1), then the usual Melnikov's theory applies. In particular 
there exists at least one subharmonic solution, and it is analytic in the perturbation parameter e. 

• If M(to) has a zero of order n, then in general no result can be given about the existence of 
subharmonic solutions. However one can introduce an infinite sequence of polynomial equations, 
which are defined iteratively: if the first equation admits a real non-zero root and all the following 
equations admit a real root, then a subharmonic solution exists, and it is a function analytic in 
suitable fractional power of e; more precisely it is analytic in. rj = e^/^, for some p < n!, and hence 
it is analytic in e^/"' . If at some step the root is simple, an algorithm can be given in order to 
construct recursively all the coefficients of the series. 

• If we further assume that the order n of the zero is odd, then we have that all the equations of 
the sequence satisfy the request made above on the roots, so that we can conclude that in such a 
case at least one subharmonic solution exists. Again, in order to really construct the solution, by 
providing an explicit recursive algorithm, we need that at a certain level of the iteration scheme a 
simple root appears. 

• Moreover we have at most n periodic solutions bifurcating from the unperturbed one with initial 
phase tQ. Of course, to count all subharmonic solutions we have also to sum over all the zeroes of 
the subharmonic Mclnikov function. 

• Finally, if M(to) vanishes identically as a function of to, then we have to extend the analysis 
up to second order, and all the cases discussed above for M{to) have to repeated for a suitable 
function Mi{to), which is obtained in the following way. If M(to) = then the solution t 
(a(t), A{t)) is defined up to first order - as it is easy to check -, so that one can expand the function 
G{a{t), A{t),t + to) up to first order: we call Mi(<o) its average over a period of the unperturbed 
solution. In particular if also Ajfi(to) vanishes identically then one can push the perturbation theory 
up to second order, and, after expanding the function G{a(t), A{t),t + to) up to second order, one 
defines M2{to) as its average over a period, and so on. 

The first conclusion we can draw is that in general we cannot say that for any vector field (F, G) there 
is at least one subharmonic solution of given period. We need some condition on G. We can require for 
G to be a zero-mean function, so that it has at least one zero of odd order. For instance, this holds true 
if the vector field is Hamiltonian, since in such a case G is the a-derivative of a suitable function. The 
same result follows if the equations describe a Hamiltonian system in the presence of small friction - how 
small depends on the particular resonance one is looking at [6,. But of course, all these conditions are 
stronger than what is really needed. 

A second conclusion is that, even when a subharmonic solution turns out to exist (and to be analytic 
in a suitable fractionary power of the perturbation parameter), a constructive algortithm to compute it 
within any given accuracy cannot be provided in general. This becomes possible only if some further 
assumption is made. So there are situations where one can obtain an existence result of the solution, but 
the solutin itself cannot be constructed. Note that such situations are highly non-generic, because they 
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arise if one finds at each iterative step a polynomial with multiple roots - which is a non-generic case; cf. 
Appendix [X] 

The methods we shall use to prove the results above will be of two different types. We shall rely on 
standard general techniques, based on the Weierstrass preparation theorem, in order to show that under 
suitable assumptions the solutions exist and to prove in this case the convergence of the series. Moreover, 
we shall use a combination of the Newton-Puiseux process and the diagrammatic techniques based on 
the tree formalism [Ul^l H] in order to provide a recursive algorithm, when possible. Note that in such 
a case the convergence of the Puiseux series follows by explicit construction of the coefficients, and an 
explicit bound of the radius of convergence is obtained through the estimates of the coefficients - on the 
contrary there is no way to provide quantitative bounds with the aforementioned abstract arguments. 
These results extend those in [5], where a special case was considered. 

The paper is organised as follows. In Section [2] we formulate rigorously the problem of subharmonic 
bifurcations for analytic ordinary differential equations in the plane, and show that, if the subharmonic 
Melnikov function admits a finite order zero, the problem can be reduced to an analytic implicit equation 
problem - analyticity will be proved in Section |4] by using the tree formalism. In Section [3] we discuss the 
Newton-Puiseux process, which will be used to iteratively attack the problem. At each iteration step one 
has to solve a polynomial equation. Thus, in the complex setting 2] the process can be pushed forward 
indefinitely, whereas in the real setting one has to impose at each step that a real root exists. If the 
order of zero of the subharmonic Melnikov function is odd, the latter condition is automatically satisfied, 
and hence the existence of at least one subharmonic solution is obtained (Theorem [1]). If at some step 
of the iteration a simple root appears, then we can give a fully constructive algorithm which allows us to 
estimate the radius of analyticity and to approximate the solution within any fixed accuracy (Theorem 
[5]) . This second result will be proved in Sections [5] and again by relying on the tree formalism; some 
more technical aspects of the proof will be dealt with in Appendix |B] Finally in Section [7] we consider 
the case in which the subharmonic Melnikov function vanishes identically, so that one has to repeat the 
analysis for suitable higher order generalisations of that function. This will lead to Theorems [3] and 21 
which generalise Theorems [5] and [TJ respectively. 



2 Set-up 

Let us consider the ordinary differential equation 

a ^ ojiA) +eF{a,A,t), 

(2.1) 

A^eG{a,A,t), 

where (a, ^) G := T x W , with C M an open set, the map A ^ oj{A) is real analytic in A, and the 
functions F, G depend analytically on their arguments and are 27r-periodic in a and t. Finally e is a real 
parameter. 

Set a{){t) = Lo[Aa)t and v4o(i) = Ag. In the extended phase-space Al x M, for e = 0, the solution 
(Q;o(i), A{){t), i-f to) describes an invariant torus, which is uniquely determined by the "energy" Aq. Hence 
the motion of the variables (a. A, t) is quasi-periodic, and reduces to a periodic motion whenever lo{Aq) 
becomes commensurate with 1. If uj{Aq) is rational we say that the torus is resonant. The parameter 
to will be called the initial phase: it fixes the initial datum on the torus. Only for some values of 
the parameter to periodic solutions lying on the torus are expected to persist under perturbation: such 
solutions are called subharmonic solutions. 

Denote by Tq{Aq) = 27r/a;(Ao) the period of the trajectories on the unperturbed torus, and define 
ui'{A) := <luj{A)/dA. If a;(Ao) = p/q G Q, call T = T{Ao) = 27rg the period of the trajectories in the 
extended phase space. We shall call p/q the order of the corresponding subharmonic solutions. 
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Hypothesis 1. One has uj'{Aq) ^ 0. 
Define 



1 '•^ 



M{to):=- dtG{ao{t),Ao,t + to), (2.2) 
^ Jo 

which is called the subharmonic Melnikov function of order q/p. Note that M{to) is 27r-periodic in to- 
Hypothesis 2. There exist to G [0,27r) and n £ N such that 

^M(io)=0 V0<fc<n-1, D(to) := ^M(to) ^ 0, (2.3) 

dig QIq 

^/lat is to is a zero of order n for the subharmonic Melnikov function. 

For notational semplicity, we shall not make explicit the dependence on t^ most of times; for instance 
we shall write D{to) ~ D. For any T-periodic function F we shall denote by (F) its average over the 
period T. 

The solution of (|2.ip with initial conditions {a{0),A{0)) can be written as 
where we have denoted by 

^'f^') (2.5) 

the Wronskian matrix solving the linearised system, and set 

^t)=eF{t)+uj{A{t))-uj{Ao)~uj'(Ao){A(t)-Ao), T{t) ^ eG{t). (2.6) 

shortening F{t) = F{a(t), A{t), t + to) and G(t) = G{a(t),A{t), t + to). 
By using explicitly (|2.5p in (|2.4p we obtain 



(2.7) 



with the notations 

In order to obtain a periodic solution we need for the mean (F) of the function F to be zero. In this 
case, if we fix also 

^^^^^ ^(^)= rdr'(F(r')-(F)), (2.8) 

^ l^oj Jo 

then the corresponding solution turns out to be periodic. So, instead of (|2.7p . we consider the system 



a{t)=a{0)+tLu'{Ao)A{0)+ f dr $(t) + c^'(^o) / dr / dr'F(T'), 
A{t)^A{0)+ /*dTF(T), 





a{t) = a{Q) +JdT {^{r) - ($}) +a;'(Ao) J dT(g(r) - (5)), 

A(0 = A{0) + g{t), ° (2-^) 

(r) = 0, 

where A{0) is determined according to (|2.8p and a(0) is considered as a free parameter. 
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We start by considering the auxiliary system 

a{t) ^ a{0) + f At (^(t) - (<&)) + ^'(Aq) / dT(g(r) - (0)), 



(2.10) 



A{t)^A{Q) + g{t), 



that is we neglect for the moment the condition that the mean of F has to be zero. Of course, only in 
that case the solution of (|2.10p is solution also of (|2.9p . hence of (|2.7p . 

It can be more convenient to work in Fourier space. As we are looking for periodic solutions of period 
T — 2'Kq, i.e. of frequency cu = 1/q, we can write 

a{t)^ao{t)+m, /3W=^e'"''%, A{t) = Aq + B{t), B{t) =Y,e^''*B,. (2.11) 
If we expand 

G(a,Ai + to) =5]5]e'''"+''''(*+*'')G,,,K^), G,y{AM) := ^'''''G^yiA), (2.12) 

with an analogous expressions for the function $(i), then we can write 

F(t) = e'^^-^'F^, $(t) = e"^""'^^, (2.13) 

with 



r- = ^EE E ^(zz^o)''aiG,„,,. (Ao, to) ,i?.„+,...B..+,, (2.14a) 

OO CXD ^ 

'^-^^EE E ^(*z^o)''aAi^.o,..^(^o,to)/3.i.../?..i?..+i...S..+. 

OO 

+ E E 79XAo)i3.,...S.,. (2.14b) 



s< 

s—2 ui + ... + 1's—i'' 



Then ((27TO| becomes 



F 



for 7^ 0, provided 



/3o = a(0) — - ^'(^o) E 7^' 



(2.15) 



(2.16) 



for 1^ = 0. Also (j2.9p can be written in the same form, with the further constraint Fq — 0. 

Then we can use /3o as a free parameter, instead of a(0). This means that we look for a value of /3o 
(depending on e) such that, by defining Bq according to the second equation in (j2.16|) . the coefficients 



6 



Bi, are given by (|2.15p for v 0. In other words, in Fourier space (|2.10p becomes 



$0 ■ ■ (2.17) 



whereas f3o is left as a free parameter. 

We look for a solution {a(t),A{t)) of (|2.10p which can be written as a formal Taylor series in e and 
(3o, so that 



oo oo 



= ait; e, M = «o(i) + E E ^''^o P (0, (2.18a) 

fe=l i=0 



oc oc 

A(t) = A{t; e, /3o) = + ^ E ^"^^o (i), (2.18b) 

fc=l j=0 

which reduces to (Q;o(i),^o) as e 0. By comparing (|2.18p with (|2.1ip we can write the Fourier 
coefficients of the solution (a(i), A(<)), for 7^ 0, as 

00 00 00 00 

:5,=^,(£,/3o)-EE-^'^o^^'''^ B. = S,(£,/3o)=EEs'/'o^^'''\ (2-19) 
k=i j=a k=i j=o 

(}^ j\ 

and, analogously, Sq ' is the contribution to order in e and j in (3q to Bq. 

By analyticity also the function r(<) = eG{a{t), A{t),t + to) can be formally expanded in powers of e 
and /3o, and one has 

00 oc 00 00 

m = T{t; £, /3o) = E E e'/J^r^"''^ (^) = E E ^'(^0 E e^-'-'f^ , (2.20) 

fc=lj=0 k=l j=0 t/GZ 

where each F^ ' is expressed in terms of the Taylor coefficients of p.l9p of order strictly less than k,j. 
By definition one has (F' = Fq The same considerations hold for $(t) = ^{a{t), A{t),t + to). 
Hence one can formally write, for all fc > 1 and j > 

P„ — — 1-^(^0)7: 7^, ^i, — — , Vi^O 

' ~ u^'iAoY 

Lemma 1. For any /3o S R the system i2.1(J\) admits a solution (a{t), A{t)) which is T-periodic in time 
and analytic in e, depending analytically on the parameter [3o- 

Proof. One can use the tree formalism introduced in Section HI see in particular Proposition [T] ■ 
It can be convenient to introduce also the Taylor coefficients 

(/3o) = E /^o fu''^ : b'^u^ (/3o) = E /^o b'^u'''' > 

(2 22) 



j>0 j>0 



Note that tI^\o) = and so on. 
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Lemma 2. Consider the system i2.10\} . Assume that To (0) = for all k eN. Then for /9o = the 
solution {a{t),A{t)) of \2.10]) is also a solution of \2.9\) . 

Proof. Simply note that ([TTOl) reduces to if T^^^ (0) = for all ken. ■ 

Of course we expect in general that rQ*^' (0) do not vanish for all A: G N. In that case, let fco G N be 
such that To''' (0) = for /c = 1, . . . , /fco and r|,''"^^^ (0) ^ 0. 
Let us define 

^(°)(s,/3o) E ^'Pi^^l -^S = (2.23) 

fe,i>0 

so that eJ^(°)(e,/3o) = (r( • ;e,/3o))- 

Lemma 3. T^'^'' (e, /3o) is Po-general of order n, i.e. JT^'^j = for j — Q, . . . ,n — I, while Tq^I ^ 0. 

Proof. This can be easily shown using the tree formalism introduced in Section [51 In fact for all j, 
■^o'j ~ ''^ associated with a tree with 1 node and j leaves. Hence one has 

jTi'^'^ = {diG{ao{-),Ao,-+to)) - (-a;(Ao))-^-^(M, (2.24) 

dig 

where the second equality is provided by Lemma 3.9 on [H]. Then .7-"*^°^ (e, /3o) is /3o-general of order n by 
Hypothesis [21 ■ 

Our aim is to find Po = (3o{e) such that JF(°^(e, /3o(e)) = 0. For such (3o a solution of (|2.10p is also 
solution of ()2.9p . If we are successful in doing so, then we have proved the existence of subharmonic 
solutions. 



3 The Newton-Puiseux process and main results 

Given a convergent power series JT^") (e, /3o) G ]R{e,/3o} as in (j2.23p . we call carrier of JF*^°^ the set 

A(^(")) := {{k,j) e N X N : .Fg ^ 0}. (3.1) 
For all V G A(jr('')) let us consider the positive quadrant 2t„ := {v} + (R+)^ moved up to v, and define 

2t:= U 2l„. (3.2) 

«eA(j^(o)) 

Let C be the convex hull of 21. The boundary dC consists of a compact polygonal path 7-"^°^ and two half 
lines Tzf^ and 71^2 ■ The polygonal path 7^^°) is called the Newton polygon of 

Notice that if the Newton polygon is a single point or, more generally, if T^^^ = for all fc > then 

there exists j > 1 such that J-"*^°^(e, /3o) = Pq ■ 5(e,/3o) with 0{e,0) ^ 0, hence /3o = is a solution of 
equation JF(°^(£, 0) = 0, that is the conclusion of Lemma [H Otherwise, if we further assume that J^(o) 
is /9o-general of some finite order n, there is at least a point of A(jr('')) on each axis, then the Newton 

polygon r^°^ is formed by iVo > 1 segments pf\ . . . and we write P'^^l = "Pf ^ U . . . U T'^^^ cf. 

Figure [H 
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Figure 1: Newton polygon. 



For all i = 1, . . . , iVo let — e Q be the slope of the segment V^^^ so that one can partition 

according to the weights given by : 

^W(.,/?o) = ^f (.,/3o) + er(£,/3o)= E ^i^f^'o + E ^S^'/^o, (3.3) 

where rl"'' is the intercept on the fc-axis of the continuation of vj;^^ . 

Hence the first approximate solutions of J^{e,Po) = are the solutions of the quasi-homogeneous 
equations 

^f)(e,/3o)= E •^S^'/5o = 0, i = l,...,No, (3.4) 

where f)i'^Vpi°^ = with ^^'^\pi°^ relatively prime integers, and s-°^ = p-^^t^*^^. 
We introduce the polynomials P/°^ = P/°^ (c) in such a way that 

^P{s,c{aoer'^') = {aosy'''' E Qkjc^ = M'''"' (c), <7o := sign (s), (3.5) 
where = ^i°]cr^- 



Lemma 4. With the notation introduced before, let 11^^ be the projection of the segment vf^ on the 
j-axis and let £i = £(11^°^) be the length ofuf^K Then P[^\c) has £i complex non-zero roots counting 
multiplicity. 



Proof. Lot m, n be respectively the maximum and the minimum among the exponents of the variable /3o 
m JC"f\ Then ii = m — n. Hence p/"' is a polynomial of degree m and minimum power n: we can write 
pI°\c) = c"P{c), where P has degree 4 and P(0) 7^ 0. Fundamental theorem of algebra guarantees that 
P(c) = has £i complex solutions counting multiplicity, which are all the non-zero roots of P/'^^ ■ 
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Let 5Ro be the set of all the non-zero real solutions of the polynomial equations P^°\c) = 0. If % = 
the system ()2.ip has no subharmonic solution, as one can easily verify. 

Let us suppose then that there exists cq G 5io, such that co(croe)^^ ' is a first approximate solution of 
the implicit equation JF^"^ (e, /3o) ~ for a suitable i — 1, . . . ,No- From now on we shall drop the label i 
to hghten the notation. We now set ei — (o-q^)^^''"" , and, as ef'"' divides JT^'^) (ctoeJ* \coef ' + yisf 
we obtain a new power series J^*-^^ (ei, yi) given by 

^(")(aoef°',co£f%yisf"') =efVW(£i,yi), (3-6) 
which is ?/i-general of order ni for some rii > 1. 

Lemma 5. With the notations introduced before, let us write P^-^\c) = go{c){c — cq)"^° with (?o(co) ^ 
and Too < n. Then ni ~ toq. 

Proof. This simply follows by the definitions of J^^^^ and P*^"^ . In fact we have 

efVW(ei,yi)=ef' [ ^ Q,,,(co + z/i^ + . .) | , (3.7) 

\/c+^(0)j=r<0) / 

SO that T<-^HO,yi) - P(")(co + yi) - <?o(co + yi)yr°, and 5o(co + yi) ^ for yi = 0. Hence ^(i) is 
j/i-general of order ni = too. ■ 

Now we restart the process just described: we construct the Newton polygon T'^^^ of T^^K If J-j^^^ ~ 
for all A: > 0, then jF(^)(ei,0) = 0, so that we have J^*^°-'(e, co((Toe)'''''') = 0, i.e. co(croe)^*''' is a solution 
of the implicit equation JF'^°) (e, /3o) = 0. Otherwise we consider the segments . . . ,P'^^ with slopes 
— l//ip'' for alH = 1, . . . , iVi, and we obtain 

where is the intercept on the fc-axis of the continuation of Hence the first approximate solutions 
of J-'-'^'^ {ei,yi) = are the solutions of the quasi-homogeneous equations 

•?i'^(£i:2/i)- E ^k,hiyi=^' t = l,...,N,, (3.9) 

where bPVpr^ — "^i^h fj-^^p^^ relatively prime integers, and s|^^ = p-"^'r-"^^ 
Thus we define the polynomials P^'^' such that 

J^^'\e^,cef) = ' ^ -^S^' = ^^'H^), (3.10) 

and we call 5Ri the set of the real roots of the polynomials P/^'. If 5Ri = 0, we stop the process as there 
is no subharmonic solution. Otherwise we caU P(i) (i.e. again we omit the label i) the polynomial which 

has a real root ci, so that cis^ is an approximate solution of the equation J-^^\si, yi) — 0. Again we 

substitute £2 = i and we obtain 

^«(£f\ci4"'+y2ef' ) = .fV(^)(e2,2/2), (3.11) 
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which is y2-general of order n2 < rii, and so on. Iterating the process we eventually obtain a sequence of 
approximate solutions 

/3o = (croe)''*"'(co + yi), yi = ^1* '(ci +^2), Z/2 = £2"(c2 + ^3), ■•• (3.12) 
where c„ is a (real) root of the polynomial P^"-* such that 

^(")(£„,C£^:'"')=£r (P(")(C) +£„(...)), (3.13) 

for all n > 0, where the functions jr(")(£„, ?/„) are defined recursively as ^^"^(e^^^l, c„e|^'^| +y„+ie|^'_^|) = 
4+i^^^"+^Hen+i, yn+i) for n > 1, with e„+i = e^'''"' and the constants r("), s("), f)("), defined 
as in the case n = in terms of a segment P^^^ of the Newton polygon of . Therefore 

^ c,iaosf°^ + c,iaosr"''+^''' + c2(aos)'^'"'+-<^'/P<")+M<^)/P<%<^' + . . . (3.14) 

is a formal expansion of (3q as a series in ascending fractional powers of doe. This iterating method is 
called the Newton-Puiseux process. Of course this does not occur if we have 3fi„ = at a certain step 
n-th, with n > 0. 

From now on we shall suppose 3fi„ ^ for all n > 0. Set also no = n. 
Lemma 6. With the notation introduced before, ifnt+i — di :— deg(P*-'-') for some i, then /i*-') is integer. 
Proof. Without loss of generality we shall prove the result for the case i ~ Q. Recall that 

.F«(0,2/i)= ^ Qfe,,(co + yi)^' =p(°Hco + 2/i), (3.15) 

fe+^(0)j=r<0) 

with r("' = ^l^^'ldQ. If do = ni, then p("' is of the form p(°)(c) = Po(c-co)"i, with Po 7^ 0. In particular 
this means that Qk.wi-i 7^ for some integer fc > with the constraint k + /i'^"^(ni — 1) = /i'"^ni. Hence 
= fc is integer. ■ 

Lemma 7. With the notations introduced before, there exists io > such that /i^*^ is integer for all i > io- 

Proof. The series are y^-general of order n^, and the and the di form a descending sequence of 
natural numbers 

n = no > do > ni > di > . . . (3.16) 

By Lemma m /i*^*^ fails to be integers only if di > n^+i, and this may happen only finitely often. Hence 
from a certain io onwards all the /i*^*-' are integers. ■ 

By the results above, we can define p := p(°^ • . . . • p*^*") such that we can write (|3.14p as 

/?o=/3o(£)= E /3o"'('^o£)''/^ (3.17) 

h>ho 

where ho — [)*^°)p(^^ • . . . • p^'"-'. By construction JT^'') (e, /3o(e)) vanishes to all orders, so that (|3.17p is a 
formal solution of the implicit equation J^^^'>{e,Po) = 0. We shall say that p.lTp is a Puiseux series for 
the plane algebroid curve defined by T^^^e.Po) = 0. 

Lemma 8. For all i > we can bound p'-'-' < n^. 
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Proof. Without loss of generality we prove the result for i = 0. By definition, there exist k',j' integers, 
with / < no, such that 

h(0) ^ , - k' 

5^ = ,» = L_A. ,3.18, 

and \)^^\lp^^^ are relatively prime integers, so that p^"^ < j' < Uq. ■ 

Note that by Lemma [S] we can bound p < no • . . . • n^j, < no!. 

Lemma 9. Let J-^^^ {e, (3o) G II^I^i/^o} be Po-general of order n and let us suppose that 5R„ ^ for all 
n > 0. Then the series \3.11^ , which formally solves J^^^^ {e, Poi^)) = 0; convergent for e small enough. 

Proof Let P^(o)(e;/3o) be the Weierstrass polynomial [2] of in C{e}[^o]- If is irreducible in 
C{£, Po}, then by Theorem 1, p. 386 in [1], we have a convergent series (3o{e^^^) which solves the equation 
:F(o)(e,/3o) = 0. Then all the following 

/3o (e^/") , Po ((e^-V)!/") , . . . , /3o ((e2-("-i)V)i/") (3.19) 

are solutions of the equation (e, /3o) = 0. Thus we have n distinct roots of the Weierstrass polynomial 
P^(o) and they are all convergent series in Cje^/"}. But also the series (|3.17l) is a solution of the equation 
T''°\e,(3o) = 0. Then, as a polynomial of degree n has exactly n (complex) roots counting multiplicity, 
(|3.17|1 is one of the p.l9p . In particular this means that (|3.17p is convergent for e small enough. 
In general, we can write 

JV 

T^°He,M^J[i:Fl"\e,f3o)r^, (3.20) 

for some iV > 1, where the JF^°'' are the irreducible factors of T^^K Then the Puiseux series (|3.17p solves 
one of the equations J-'^^\e, Po) = 0, and hence, by what said above, it converges for e small enough. ■ 

As a consequence of Lemma [9] we obtain the following corollary. 

Theorem 1. Consider a periodic solution with frequency lo = p/q for the system Assume that 

Hypotheses [7] and [H are satisfied with n odd. Then for e small enough the system \2.1\) has at least one 
subharmonic solution of order q/p. Such a solution admits a convergent power series in \e\^^"' , and hence 
a convergent Puiseux series in \e\. 

Proof. If n is odd, then 3fi„ 7^ for all n > 0. This trivially follows from the fact that if n is odd, then 
there exists at least one polynomial P^^'' associated with a segment T'l"'' whose projection Tlf^"^ on the 

j-axis is associated with a polynomial Pj^^^ with odd degree £i . Thus such a polynomial admits a non-zero 
real root with odd multiplicity ni, so that T'-^\ei,yi) is j/i-general of odd order ni and so on. 

Hence we can apply the Newton-Puiseux process to obtain a subharmonic solution as a Puiseux series 
in £ which is convergent for e sufficiently small by Lemma [HI ■ 

Theorem [1] extends the results of [T3] . First it gives the explicit dependence of the parameter Pq on 
£, showing that it is analytic in |£|^/"'. Second, it shows that it is possible to express the subharmonic 
solution as a convergent fractional power series in e, and this allows us to push perturbation theory to 
arbitrarily high order. 

On the other hand, the Newton-Puiseux algorithm does not allow to construct the solution within 
any fixed accuracy. In this regard, it is not really constructive: we know that the solution is analytic in 
a fractional power of £, but we know neither the size of the radius of convergence nor the precision with 
which the solution is approximated if we stop the Newton-Puiseux at a given step. Moreover we know 
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that there is at least one subharmonic solution, but we are not able to decide how many of them are 
possible. In fact, a subharmonic solution can be constructed for any non-zero real root of each odd-degree 
polynomial P/"'' associated with each segment of P*^") to all step of iteration, but we cannot predict a 
priori how many possibilities will arise along the process. 

However, we obtain a fully constructive algorithm if we make some further hypothesis. 

Hypothesis 3. There exists Jq > such that at the iQ-th step of the iteration, there exists a polynomial 
pM — pM(^f,'j mJiiQfi jidg a simple root Ci„ £ R. 

Indeed, if we assume Hypothesis [31 we obtain the following result. 

Theorem 2. Consider a periodic solution with frequency oj = p/ q for the system h2.1\l . Assume that 
Hypotheses\^\^and\^are satisfied. Then there exists an explicitly computable value > such that for 
\e\ < Eq the system k2.1\) has at least one subharmonic solution of order q/p. Such a solution admits a 
convergent power series in jej"'^/"', and hence a convergent Puiseux series in 

We shall see in Section [5] that, by assuming Hypothesis [31 we can use the Newton-Puiseux algorithm 
up to the ig-th step (hence a finite number of times), and we can provide recursive formulae for the 
higher order contributions. This will allow us - as we shall see in Section [H] - to introduce a graphical 
representation for the subharmonic solution, and, eventually, to obtain an explicit bound on the radius 
of convergence of the power series expansion. 

4 Trees expansion and proof of Lemma [1] 

A tree 9 is defined as a partially ordered set of points [vertices) connected by oriented lines i. The 
lines are consistently oriented toward a unique point called the root which admits only one entering line 
called the root line. If a line i connects two vertices 0i,D2 and is oriented from 02 to Oi, we say that 
02 -< Oi and we shall write = ^- We shall say that i exits 02 and enters Di. More generally we write 
X>2 -< Oi when Oi is on the path of lines connecting O2 to the root: hence the orientation of the lines is 
opposite to the partial ordering relation -<. 

We denote with V{0) and L(6') the set of vertices and lines in 9 respectively, and with |V^(0)| and 
|L(0)| the number of vertices and lines respectively. Remark that one has = \L(d)\. 

We consider two kinds of vertices: nodes and leaves. The leaves can only be end-points, i.e. points 
with no lines entering them, while the nodes can be either end-points or not. We shall not consider the 
tree consisting of only one leaf and the line exiting it, i.e. a tree must have at least the node which the 
root line exits. 

We shall denote with N{9) and E{d) the set of nodes and leaves respectively. Here and henceforth we 
shall denote with and e the nodes and the leaves respectively. Remark that V{d) ~ N{9) H E{6). 

With each line £ = io, we associate three labels {h£, 6i, vi), with h( e {a. A], 5( G {1, 2} and e Z, 
with the constraint that vg ^ for hi = a and Sg = 1 for hi — A. With each line £ = we associate 
hi = a, 5i ^ I and i/f = 0. We shall say that hg, 6g and vg are the component label, the degree label and 
the momentum of the line £, respectively. 

Given a node 0, we call r„ the number of the lines entering carrying a component label h — a and 
St, the number of the lines entering with component label h — A. We also introduce a badge label 
60 S {0, 1} with the constraint that &o = 1 for hg^ = a and 5g^ — 2, and for hg^ — A and ug^ 7^ 0, and 
two mode labels ctdj S 2- We call global mode label the sum 

J^D = pern + gcTo, (4.1) 
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where q,p are the relatively prime integers such that oj{Aq) = p/q, with the constraint that Vo = when 
K=0. 

For all £ = If,, we set also the following conservation law 



(4.2) 



to e N(0) 
ro -< D 



i.e. the momentum of the line exiting is the sum of the momenta of the lines entering plus the global 
mode of the node itself. 

Given a labeled tree 6, where labels are defined as above, we associate with each line £ exiting a node, 
a propagator 



{iuJi^eY' 
1 



while for each line £ exiting a leaf we set ge = 1. 

Moreover, we associate with each node a node factor 



, he = a, A, ve ^ 0, 
hi = A, Vi = 0, 



(4.3) 



K = < 



rt,\s„\ 

^u;{Ao), 

ro!so! 



ro!so! 



32 



he. 


= a, 




= 1, 


K 


= 1, 


i^e, 7^ 0, 


he. 


= a, 




= 1, 


K 


= 0, 


J^e, + 0, 


he. 


= OL, 


k. 


= 2, 


K 


= 1, 


n, 0, 


he. 


= A 


Se, 


= 1, 


K 


= 1, 


vi, 0, 


he. 


= A 


Se, 


= 1, 


K 


= 1, 


yi, = 0, 


he. 


= A 




= 1, 


K 


= 0, 


yi, = 0, 



(4.4) 



with the constraint that when feo = one has rt, = and So > 2. 

Given a labeled tree Q with propagators and node factors associated as above, we define the value of 
Q the number 



Val(^) = I n 5^ 

ye(kL{6) 



iDeiv(6i) 



(4.5) 



Remark that Val(^) is a well-defined quantity because all the propagators and node factors are bounded 
quantities. 

For each line £ exiting a node we set he = h^, while for each line £ exiting a leaf we set 6^ = 0. Given 
a labeled tree 0, we call order of Q the number 



k{Q) = \{£^L{Q) ■ 6, = 1}|; 



(4.6) 



the momentum v{S) of the root line will be the total momentum, and the component label h(&) associated 
to the root line will be the total component label. Moreover, we set j{9) = \E{9)\. 

Define %^u,h,j the set of all the trees 9 with order k{9) = k, total momentum i'{9) = v, total component 
label h{9) = h and j{9) = j leaves. 
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Lemma 10. For any tree 9 labeled as before, one has \L{9)\ = < 2k{d) + j{0) — 1. 

Proof. We prove the bound < 2k{d) — 1 by induction on k. 

For k — 1 the bound is triviaUy satisfied, as a direct check shows: in particular, a tree 9 with k{9) = 1 
has exactly one node and j{d) leaves. In fact if 6 has a line i = ia with bg = 0, then has s„ > 2 lines 
with component label h ^ A entering it. Hence there are at least two lines exiting a node with b„ = 1. 

Assume now that the bound holds for all k' < k, and let us show that then it holds also for k. Let £o 
be the root line of 9 and Oq the node which the root line exits. Call r and s the number of lines entering 
Oq with component labels a and A respectively, and denote with 6*1 , ... , 9r+s the subtrees which have 
those lines as root lines. Then 

r+s 

\N{9)\ = l+Y,mO„.)l (4.7) 

m— 1 

If £o has badge label bg^ = 1 we have \N{9)\ < 1 + 2{k - I) - {r + s) < 2k - 1, by the inductive 
hypothesis and by the fact that k{9i) + . . . + k{9r+s) — k — 1. li £o has badge label big = we have 
W{9)\ <l + 2k-{r + s) < 2fc- 1, by the inductive hypothesis, by the fact that k{9i) + . .. + k{9r+s) = k, 
and the constraint that s > 2. Therefore the assertion is proved. ■ 

Lemma 11. The Fourier coefficients [3^^' , v ^ Q, and ' can be written in terms of trees as 



(3l-'>= Val(^)' '^^O' (4.8a) 

-^ik.j) ^ j2 Val(0), 1/ e Z, (4.8b) 



for all k>l,j>0. 

Proof. First we consider trees without leaves, i.e. the coefficients /J^ ' , ly ^ 0, and ' . For fc = 1 is a 
direct check. Now let us suppose that the assertion holds for all k < k. Let us write fa = P, fA = B and 
represent the coefficients f^'^if^ with the graph elements in Figure [H as a line with label v and h = a, A 
respectively, exiting a ball with label (fc,0). 



(fc,0) (fc,0) 



(fc,0) 



A 1 



Figure 2: Graph elements. 

Then we can represent each equation of (|2.2ip graphically as in Figure [31 simply representing each 

(k- 0) r?n 

factor ff^_''l. in the r.h.s. as a graph element according to Figure[21 the lines of all such graph elements 
enter the same node Oq. 
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h s 



- r 




Figure 3: Graphical representation for tlie recursive equations 1 12. 231 1. 



The root line £o of such trees will carry a component label h — a,A for f ^ P,B respectively, and a 
momentum label v. Hence, by inductive hypothesis, one obtains 



Jh,v 



2^ 9io-^^ofh^ 



(fci,0) 



E Val(0) 

,1^1 , hi ,0 



E 



Val(6l) 



E Val(^), 



(4.9) 



where m = + sq, and we write ^* for the sum over all the labels admitted by the constraints, so that 
the assertion is proved for all k and for j ~ 0. 

Now we consider k as fixed and we prove the statement by induction on j. The case j = has already 
been discussed. Finally we assume that the assertion holds for j — j' and show that then it holds for 
j' + 1. Notice that a tree 6 G Tk^u.h,j'+i, for both h = a,A can be obtained by considering a suitable 
tree Oq £ Tic^u,h,j' attaching an extra leaf to a node of 6q and applying an extra derivative da to the node 
factor associated with that node. If one considers all the trees that can be obtained in such a way from 
the same a-nd sums together all those contributions, one finds a quantity proportional to daVal{9o)- 

Then if we sum over all possible choices of Oq, we reconstruct (3^ ' for h — a and ' for h = A. 
Hence the assertion follows. ■ 



Proposition 1. The formal solution \2.19\) of the system i2.10\) . given by the recursive equations h2.21]) . 
converges for e and /3o small enough. 

Proof. First of all we remark that by Lemma 1101 the number of unlabeled trees of order k and j leaves 
is bounded by 4^*^+^ x 2^'^+-' — 8^*^+^ . The sum over all labels except the mode labels and the momenta 
is bounded again by a constant to the power k times a constant to the power j, simply because all such 
labels can assume only a finite number of values. Now by the analyticity assumption on the functions F 
and G, we have the bound 



(4.10) 
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for suitable positive constants Q, R, S, k, and we can imagine, without loss of generality, that Q and S 
are such that \d\uj{Ao)/ s\\ < QS'^. This gives us a bound for the node factors. The propagators can be 
bounded by 

\gi\ < max 

so that the product over all the lines can be bounded again by a constant to the power k times a constant 
to the power j. 

Thus the sum over the mode labels - which uniquely determine the momenta - can be performed by 
using for each node half the exponential decay factor provided by (|4.10p . Then we obtain 

< CiC2^C|e-«l''l/2, \Bi^''^\ < ClC2*C^e-'^l'^l/^ (4.12) 

for suitable constants Ci, C2 and C3. This provides the convergence of the series (|2.19p for |e| < Cj"^ 
and 1^0 1 < ^3^^ ■ 



cj'iAo) 



1 
to 



(4.11) 



5 Formal solubility of the equations of motion 

Assume that Hypotheses [H [5]and[3]are satisfied. Let us set rj :— \e\^^'^, where p = p^''^ • . . . • p^'^°\ We 
search for a formal solution {a{t),A{t)) of ([2T|) . with a{t) = ao{t) + /Sq + P{t) and A(<) = A^ + B{t), 
where ^ ^ 

fe>i i^ez fe>i vei k>i 

and the coefBcients /3g'^' , /Si*^' and bI^^ solve 



^[k] Uk] 

lUJl/ (ILUV)'^ IUJI> 

[k] 



_ __5o_ rW-n 
with the functions T^!f^ and $1*^' recursively defined as 



(5.2) 



r! s. 

m>0 r+s=m pcro+(3(To+i/i + ...+!^„=!^ 
fci+...+fcm=A:— p 

^l'' = E E E ^f^.o.i(^o,to)/3[^l•••/3[^l4^-'•••4^^ 

/ci + ... + A;m— A;— p 

+ E E f -(^0)4^1... (5.3b) 

s>2 i^i + ... + i/s = i/ 

/£l + ... + fea=fc 

where /Slf' = /3?1 for 1/ 7^ . We use a different notation for the Taylor coefficients to stress that we are 
expanding in 77. 

We say that the integral equations ()2.9p . and hence the equations (|5.2p . are satisfied up to order k if 

there exists a choice of the parameters /3q^' , . . . /Jq*^' which make the relations (|5.2p to be satisfied for all 
k 1 , . . . , . 
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Lemma 12. The equations 15. are satisfied up to order k = p — 1 with (31 and Bl identically zero 
for all k — 1, . . . , p — 1 and for any choice of the constants /3q^' , . . . (3^^ . 

Proof. One has e ~ <J'if , with a = sign(e), so that $1''' = r[f'' = for aU fc < p and all G Z, 
independently of the values of the constants (3^q \ . . . (3^^ Moreover /^if' = slf' = for all fc < p. ■ 

Lemma 13. The equations / f 5. 2fl are satisfied up to order k — p, for any choice of the constants 

Proof. One has = G(ao(t), ^o, * + ^o) and $[Pl = F(ao(i), Aq, i + to), so that 

4"!= E G<,„...(Ao), $1^1= E F'^o.^^Ao). (5.4) 

Thus, /^if' and si""' can be obtained from jO]). Finally r|f' = Af(io) by definition, and one has Af (io) = 
by Hypothesis [51 Hence also the last equation of (|5.2p is satisfied. ■ 

Let us set 



f)o = 


(,(0) 




. . . • p*-*"^ 




So 


= s(o) 


p(l'....-p(*°) 


1 


f)l = 






|p(2).... 


• p(*«). 


Si 


= So - 


hs(i)p(2)..... 


p(«o) 




t)l- 




'p(3).... 


• p(*«). 


S2 


= Si - 


t-s(2)p(3)..... 


p(*o) 


f)»o = 




-1 + 






Sio 


= Sig 


_i H-s^*"). 





(5.5) 



Lemma 14. T/ie equations i5.2\) are satisfied up to order k = p + Sig provided /3q = Ci, with Ci the 
real root of a polynomial p(*)(c) of the i-th step of iteration step of the Newton-Puiseux process, for 
i — 0, . . . ,io, and (3^^ ' = for all fc' < \)ig, fc' ^ \)i for any i. 

Proof. If /3|,'''l = for all 1< fc' < [)o, one has V^'o^ = for aU p < fc < p + 5o, while r[f+"'l = P^^\(3^^°^), 
so that r{}'+'"l = for = co. Thus r[fl = for p + Sq < fc < p + Si provided (3^^'^ = for aU 

[)o < fc' < t)i, while 4"+'^! = P(i)(/3[''^1), so that 4+'^' for Z?'"^' = ci, and so on. 

Hence if we set /3g'''' = for alH = 0, . . . , io, and /jlf ' = for all fc' < , fc' 7^ [)i for any i = 0, . . . , io, 
one has Fg — for all p < fc < p + . Moreover, and Fi;, are well-defined for such values of fc. 
Hence (|5.2p can be solved up to order fc = p + , indipendently of the values of the constants (3^^ ' for 
fc'>[}.o- ■ 

By Lemma [T4l we can write 

/3o = /3o(?7) COT?"" + cir?"^ + . . . + Q^Ty^'o + ^ /jl^'o+'^'lyyl-o+fc. (5.6) 

fc>i 

Lemma 15. The equations 115. 2\) are satisfied up to any order k = p + Si„ + k, k > 1 provided the 
constants (3^^"'°^'" ' are suitably fixed up to order k' = k. 
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Proof. By substituting (|5.6p and e = arf , with a = sign(e), in ro(e,/3o) we obtain 
^oicrv'^ , Mv)) ^ crr]'^ ^ QsujV""^ X! J{j,mo, . . . ,m^„,m) x 



si,j>0 mo + ...+mig+m=j 

m^rrii >0 



X rj 



where Qsi,j — •^iij^'*^ ^-nd 



n>0 ni + ...+n„i=n 

ni>l 



For any k > 1 one has, by rearranging the sums, 



m,mi,n,si,j>0 ni + ...+n,„=K,~n 

sip+mot)a + ...+migt)ig+mt)ig=Sig+n 

SO that the last equation of ()5.2p gives for k > 1 



mi,sij'>0 
mo + ...+miQ+l=j 

sip+mof)o + ---+"iio flio=5io 



+ 5: jam„,...,rn.„,HQ.,c-...c::- 

iTii,si,j>0,m>2 ni+...+n,„=K (5.10) 

mo + .-.-i-mi^+m— j l<ni</i—l 
sip+mof)o + ---+'nijj l)iQ+mf)ij,=5iQ 

+ ^ J(j,mo,...,m.„,m)g,„,cr...cr:- ^ . . . = 0, 

m,mi,si,j>0,n>l ni+...+nm = K,—n 

mo + ...+Tnig+m=j ni>l 

sip+mol)o + ---+"iio flio+™')io=^»o+" 

where ah terms but those in the first hue contain only coefficients /3q'''°^'' ' with k' < k. 
Recall that by Hypothesis [3] 

2_] m,„J(j,mo,...,mio,TO)Qsijc"V..c^LV'c^''' = — ^— (c,J =: C ^ 0, (5.11) 

mo + ---+"iio=J 

Sip+ITlol)0 + ---+niio l)io=5io 

so that we can use (jS.lOp to express p]^^"'^'^^ in terms of the coefficients /3q'''"~'''^ ' of lower orders k' < k. 

Thus we can conclude that the equations (|5.2p are satisfied up to order k provided the coefficients /3q'''°^'' ' 
are fixed as 

/S^"^"'' = -^G[«'l(co, . . . , c.„, /3^°+^ . . . , /3^°+'''^^l), (5.12) 

for all 1<k' <K, where GM(co, • • ■ , c^g, . . .,(3^^'°^'"'^^) is given by the sum of the second and 

third lines in ■ 

We can summarise the results above into the following statement. 

Proposition 2. The equations \5/0jl are satisfied to any order k provided the constants pI^^ are suitably 
fixed. In particular = slf' = Sq*^' = for k < p and pI^^ — for k < \)ig , k ^ \)i for any i = 0, . . . , zq- 



(5.7) 



J(j,mo, . . . ,m,,o,TO) := ; '■ — (5.8) 



^^[p+.,,+«]^ E J(j,-o,...,m.o,HQ.,,cr...c"- E/3^°+-l.../3t^-+""l, (5.9) 
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6 Diagrammatic rules for the Puiseux series 



In order to give a graphical representation of the coefficients /Jq*^' , /3if' and ' in (|5.ip , we shall consider 
a different tree expansion with respect to that of Section [H We shall perform an iterative construction, 
similar to the one performed through the proof of Lemma [Til starting from equations l|5.2p for the 
coefficients bH^^ for fc > p, and from (l5l^ for k>\),„+ 1. 

Let us consider a tree with leaves. We associate with each leaf e a leaf label — 0, ■ ■ ■ ,io. 

For fc = p we represent the coefficients /3|f' and i?if' as a line exiting a node, while for fc = 
i = 0, . . . , io we represent /Sg'''' as a line exiting a leaf with leaf label Oi. 

Now we represent each coefficient as a graph element according to Figure [H as a line exiting a ball 
with order label fc, with fc > tji^ + 1 for the coefficients /3q , and fc > p + 1 for the coefficients (31 and 
-BI^'; we associate with the line a component label hi € {f3Q,P,B}, a degree label Si £ {1,2} with the 
constraint that Si — 1 for hi — B,(3o, and momentum label i>i G Z, with the constraint that vi for 
hi = P, while vi = for /i^ = /3o- 

= . Q [k] 

00 1 ^ 



[fc] 



[fc] 



P 1 



/3 2 1/ 



bI'^ = 



B 1 



[fc] 



Figure 4: Graph elements. 

Hence we can represent the first three equations in (j5.2p graphically, representing each factor /3|>^'' 
and si'^'' in (|5.3p as graph elements: again the lines of such graph elements enter the same node Vq. 

We associate with Oq a badge label 6t,(, G {0, 1} by setting bc^ = 1 for /i^^ = (3 and ^^(, = 2, and for 
/i^Q = B and i/^^ 7^ 0. We call rtn, the number of the lines entering Oq with component label h = /3o,/3, 
and Soo the number of the lines entering Oq with component label h = B, with the constraint that if 

= one has r„g = and So„ > 2. Finally we associate with Oq two mode labels a^g^ a^'^ G Z and the 
global mode label Vx^o defined as in (|4.ip . and we impose the conservation law 

i=l 

where ti, . . . ,lr„ +s„ are the lines entering Oq. 
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We also force the following conditions on the order labels 

ki — k — p, bag = 1, 



Do ^-^f 

E 

1=1 



(6-2) 

^ ^ ki — fc, 6do ~ Oj 



i=l 



which reflect the condition on the sums in 

Finally we associate with = Oq a node factor JV^ = a^'Af^,, with a = sign (e) and J\f„ defined as in 
(|4.4p . and with the line £ = £„g a. propagator = g^, with g^ defined as in (|4.3p . The only difference 
with respect to Section [His that the component label can assume the values /3, B, which have the role of 
a, A respectively. 

The coefficients /3q'''° k> 1, have to be treated in a different way. 

j \ _____ 

First of all we point out that also the coefficients Fq ' in (|2.20p can be represented in terms of sum 
of trees with leaves as in Section |H In fact we can repeat the iterative construction of Lemma fTT| simply 
by defining j 
6oo = 1 and 



(fc j'^ 

by defining Tk^^rj as the set of the trees contributing to Fg ' , setting gi^ — 1, hig — T, Si^ — 1, vi^ — 



and no further difficulties arise. 



= l!^^£)p!^G.„„ (6.3) 

(0) _-p(si+lj). 



Recall that the coefficients Qsi.j in (|5.10p are defined as Qs,,j = -^.ij^'' = f^"' so that 

Qs^,j='^'' E Val(^). (6.4) 

Hence the summands in the second and third lines in (|5.10p can be imagined as "some" of the trees in 
T^j+i.o.Tj where "some" leaves are substituted by graph elements with hg = Pq. More precisely we shall 
consider only trees 9 of the form depicted in Figure [5l with si + 1 nodes, sq = so,o + . . . + So,io leaves, 
where so,a is the number of the leaves with leaf label a, and Sg graph elements with hg — (3q, such that 

Sip + ^ so,i:f)i + s'q^io = St„ + n, 
. (6.5) 



E 



h = (sq - l)t)jo + k-n, 



for a suitable < rt < fc — f}^,-, , with the constraint that when n ~ one has Sg > 2. We shall call £i the 
Sg lines with hg. = f3o. Such conditions express the condition on the sums in the second and third lines 
in ((5lU)l . 

The propagators of the lines exiting any among the si + 1 nodes and the node factors of the nodes 
(except the root line and the node which the root line exits) are g^ = gi and A/'„* — cr^'Mo with the 
component labels assuming the values /3, B, which have the role of a, A, respectively. We associate with 
the root line a propagator 

9l = (6.6) 

where C is defined in (|5.1ip . while the node = Oq which the root line exits will have a node factor 
= J\ft,a as in (|6.3p Finally we associate with each leaf e a leaf factor J\f* = c„,. 
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We now iterate such a process until only nodes or leaves appear. We shall call allowed trees all the 
trees obtained in such a recursive way, and we shall denote with QkM,h the set of allowed trees with order 
k, total momentum v and total component label h. 

Given an allowed tree 6 we denote with N{6), L{9) and E{6) the set of nodes, lines and leaves of 
9 respectively, and we denote with Ea{0) the set of leaves in 9 with leaf label a. We point out that 
E[9) = £;o(6') n . . . n E^^{9). We shall define the value of 9 as 

var(0)= I n A f n K] f n ■ (^.7) 

YeL{e) I \Dew(0) / \ce£;(0) / 

Finally, we denote with A(6') the set of the lines (exiting a node) in 9 with component label h — /3o 
and with N*{9) the nodes with b„ = 1; then we associate with each node in N*{9), with each leaf in 
Ea{9) and with each line in A{9) a weight p, \)a and Iji^ — p — Sig, respectively, and we call order of 9 the 
number 

k{9) = p\N*{9)\ + ((,,„ - p - s,„)|A(0)| + ^ [)c|i?a(e)|. (6.8) 

a=0 

Note that - p — Si„ < 0. 

Lemma 16. The Fourier coefficients pl^^ and can be written in terms of trees as 

Po^= E Var(^), A;>f],„ + 1, (6.9a) 

3^= J2 Var(^), fc>p, (6.9b) 

B W = Var(6l), k>p. (6.9c) 

\k] ^\k] 

Proof. We only have to prove that an allowed tree contributing to the Fourier coefficients /?g , (31 
and bI^^ has order k. We shall perform the proof by induction on k > t)ig + I for the coefficients /Jg*^', 

— \k] \k] — [pi 

and A: > p for (31 and Bl . Let us set fj^ — (3 and fs = B. An allowed tree 9 contributing to fj^l 

has only one node so that k{9) = p, while an allowed tree 9 contributing to /3q''*^^', i = 0, . . . ,«o has 
si + 1 nodes, sq = sq.o + . . . + so,io leaves, and one line A{9), and via the conditions (|6.5p we have 
Sip + so,of]o + . • . + so,iof]io = + 1' so that fc(0) = i)ig + 1. 
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Let us suppose first that for all k' < k, an allowed tree 6' contributing to ' has order k{9') = k' . 
By the inductive hypothesis, the order of a tree 9 contributing to /Jq'^' is (we refer again to Figure [5] for 
notations) 

k{e) = (si + l)p + ^so,,[}, + ^fc, + f)^o - p-5,,, (6.10) 

1=0 i=l 

and via the conditions in (j6.5p we obtain k{6) = k. 

Let us suppose now that the inductive hypothesis holds for all trees 6' contributing to fjl^'^ , k' < k. 

An allowed tree contributing to /j^*^]^ is of the form depicted in Figure where so.a is the number of 
the lines exiting a leaf with leaf label a and entering Oq, si is the number of the lines exiting a node and 
entering Oq, and s'q,s'i are the graph elements entering Oq with component label /3o and either (3 or B, 
respectively. 

So 




Sl 



Figure 6: An allowed tree contributing to f\^\^- 

If &t,„ = 1, by the inductive hypothesis the order of such a tree is given by 

k{e) = (si + l)p + ^ so + (6.11) 

1 = 4=1 

and by the first condition in (|6.2p we have k(9) — k. Otherwise if b„g = 0, we have sq + Sq — 'And, by 
the inductive hypothesis, 

k{0) = sip + J2h^k, (6.12) 

4=1 

via the second condition in (j6.2p . ■ 
Lemma 17. Let q := min{[)o,p} and let us define 

M = 2^ + 3. (6.13) 

q 

Then for all 9 G Qk,u,h one has 

\L{9)\ < Mk. (6.14) 

As the proof is rather technical we shall perform it in Appendix |B] 

The convergence of the series (|5.1|) for small 77 follows from the following result. 
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Proposition 3. The formal solution h5.1\) of the system 112. 9\) . given by the recursive equations \5.'2\) and 
I15.12\) . converges for rj small enough. 

Proof. By Lemma [T71 the number of unlabeled trees of order k is bounded by 4*^*^. Thus, the sum over 
all labels except the mode labels and the momenta is bounded by a constant to the power k because all 
such labels can assume only a finite number of values. The bound for each node factor is the same as in 
Proposition [TJ while the propagators can be bounded by 

\gi I 5: max 

so that the product over all the lines can be bounded again by a constant to the power k. The product 
over the leaves factors is again bounded by a constant to the power k, while the sum over the mode labels 
which uniquely determine the momenta can be performed by using for each node half the exponential 
decay factor provided by ()4.10p . Thus we obtain 

IPl^^ < CiC|e-"l"l/2, < CiC|e-«l"l/2, (6.16) 

for suitable constants Ci and C2. Hence we obtain the convergence for the series ()5.ip . for \r]\ < €2^. ■ 
The discussion above ends the proof of Theorem O 



1 



1 

C 



(6.15) 



7 Higher order subharmonic Melnikov functions 



Now we shall see how to extend the results above when the Melnikov function vanishes identically. 

We are searching for a solution of the form {a{t),A{t)) with a{t) = ao{t) + /3o + P{t) and A{t) — 
Aq + B{t), where 

^W-E^""*^-(^'^o)' Bit)=J2e''"^'B,i^^Po). (7.1) 



First of all, we notice that we can formally write the equations of motion as 

Pu (Po) = : hw(Ao)— (Po) - 



lOJV 



4\m 



(7.2) 



(7.3) 



Lu'iAo) ■ 

where the notations in (|2.22p have been used, up to any order fc, provided 

ro(e,/3o) =0. 

If M(to) vanishes identically, by ((T^ we have Tq^'^^ = for aU j > 0, that is T^'ilSo) = 0, for aU 
Po, and hence To{e,Po) = s'^J^^'^'>{e,Po), with ^^^^ a suitable function analytic in e, f3o. 

Thus, we can solve the equations of motion up to the first order in e, and the parameter Pq is left 
undetermined. More precisely we obtain 



= £/3^''+e/3(i)(e,/3o), 



sB 



■(1) 



eBi^\e,f3o) 



(7.4) 



5(1) ^(1) 



where B]^"' solve the equation of motion up to the first order in e, while l3i^\ b'^}^ are the corrections 
to be determined. 
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Now, let us set 

Mo(to) = M{to), Afi(io) = (0,<o), (7.5) 

where F^''^ (/3o, to) = ^''u\Po), we are stressing the dependence of F^,*^'"'^ on to- We refer to Mi (to) as 
the second order subharmonic Melnikov function. Notice that Mo(to) = rQ^''(0, to)- 

If there exist to G [0, 27r) and ni G N such that to is a zero of order ni for the second order subharmonic 
Melnikov function, that is 

^Mi(to) = VO<A:<ni-l, D = D{to) := ;^Afi(to) ^ 0, (7.6) 
dtp aiQ 

then we can repeat the analysis of the previous Sections to obtain the existence of a subharmonic solution. 
In fact, we have 

T(^\e,Po) ■■= E ^'^o^Sh - rr''^'(^o), (7.7) 

k,j>0 

where tg has to be fixed as the zero of Mi (to), so that, as 

(-..(Ao))-^-^Mi(to) = j!r|,''-''^(to), (7.8) 
dto 

for all j, as proved in [6| with a different notation, we can construct the Newton polygon of J-^'^\ which 
is /3o-general of order Ui by (|7.6p . to obtain B^^^ and f3o as Puiseux series in e, provided at each step 
of the iteration of the Newton-Puiseux algorithm one has a real root. 

(2) 

Otherwise, if Mi(to) vanishes identically, we have Fq (/3o) = for all f3o, so that we can solve the 
equations of motion up to the second order in e and the parameter /3o is still undetermined. Hence we 

set M2(to) = Fq (0,to) and so on. 

In general if Mfc'(to) = 0, for all fc' = 0, . . . , k — 1, we have Fo(e,/3o) = e*^ ^^'^ H^jPo), so that we can 
solve the equations of motion up to the K-th order in e, and obtain 

(3, = + . . . + e^pi"'' + e'^M'^ (£, /3o), (7.9a) 
= sB^^^ + ... + e-B^;^ + e^Bi'^He, /3o), (7.9b) 

where b''^ \ k' = 0, . . . , At — 1 solve the equation of motion up to the K-th order in e, while 

are the correction to be determined. 
Hence we can weaken Hypotheses [U and [3] as follows. 

Hypothesis 4. There exists k > such that for all k' — 0, . . . , k — 1, Mk'{to) vanishes identically, and 
there exist to £ [0, 27r) and n £ N such that 

d" 

— M«(to)-0 V0<j<n-1, D^D{to):^—AUto)^0, (7.10) 
dt^o dtS 

that is to is a zero of order n for the K-th order subharmonic Melnikov function. 

Hypothesis 5. There exists io > such that at the io-th step of the iteration of the Newton-Puiseux 
algorithm for there exists a polynomial p('") = p(*°)(c) which has a simple root c* £ M. 

Thus we have the following result. 
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Theorem 3. Consider a periodic solution with frequency uj = p/q for the system \2.1\) . and assume that 
Hypotheses\^^and\^are satisfied. Then there exists an explicitly computable value Eq > such that for 
\e\ < Eo the system h2.1\) has at least one subharmonic solution of order q/p. Such a solution admits a 
convergent power series in \e\^^"' , and hence a convergent Puiseux series in \e\. 

The proof can be easily obtained suitably modifying the proof of Theorem [21 

Now, call ^it^ the set of real roots of the polynomials obtained at the n-th step of iteration of the 
Newton-Puiseux process for J'^^"). A gain if n is even we can not say a priori whether a formal solution 
exists at all. However, if SRI'"'' ^ for all rt > 0, then we obtain a convergent Puiseux series as in Section 

13 

Finally, as a corollary, we have the following result. 

Theorem 4. Consider a periodic solution with frequency u = p/q for the system !i2.1]) . Assume that 
Hypotheses [7] and are satisfied with n odd. Then for e small enough the system 12.1]) has at least one 
subharmonic solution of order q/p. Such a solution admits a convergent power series in |£|"', and hence 
a convergent Puiseux series in e. 

Again the proof is a suitable modification of the proof of Theorem [TJ 
Acknowledgements. We thank Edoardo Sernesi for useful discussions. 



A On the genericity of Hypothesis 3 

Here we want to show that Hypothesis [3] is generic on the space of the coefficients of the polynomials. 
More precisely, we shall show that given a polynomial of the form 



Pia,c) = Y, 



n > 1, 



a := (ao, . . . ,a„), 



(A.l) 



i=0 



the set of parameters (ao,...,a„) G M"+^ for which P{a,c) has multiples roots, is a proper Zariski- 
closecQ subset of K"+^. Notice that a polynomial P = -P(a, c) has a multiple root c* if and only if also 
the derivative dP/dc vanishes at c* . 
Recall that, given two polynomials 



Pl{c) = X! P2(c) = ^ b,n~iC\ 

i=0 i=0 



(A.2) 



with n,m > 1, the Sylvester matrix of Pi, P2 is an n + to square matrix where the columns 1 to m are 
formed by "shifted sequences" of the coefficients of Pi , while the columns m + 1 to to + n are formed by 
"shifted sequences" of the coefficients of P2, i-e. 





. 





bo 


. 


■ \ 


a-i 


ao . 





hi 


60 . 








. 







. 


■ bm~l 




. 


On 





. 


bm J 



Syl(Pi,P2) := 



and the resultant P(Pi,P2) of Pi,P2 is defined as the determinant of the Sylvester matrix. 



(A.3) 



^See for instance 1121. 



26 



Lemma 18. Let ci.i, . . . , ci^n o.'^d C2,i, ■ . . , C2.m the complex roots of Pi, P2 respectively. Then 

n rn 

i?(Pi,P2) =ao"'&onn(^M-C2,,). (A.4) 

A complete proof is performed for instance in [13j . In particular, Lemma [T51 implies that two polyno- 
mials have a common root if and only if R{Pi, P2) =0. 

Recall also that given a polynomial P — P{c), the discriminant D{P) of P is the resultant of P and 
its first derivative with respect to c, i.e. D{P) := R{P,P'), where P' :— dP/dc. Thus, a polynomial 
P — P{a, c) of the form (jA.ip has a multiple root if and only if its discriminant is equal to zero. 

Now let us consider the set 

F {a = (flo, . . . , a„) e M"+^ : P{a, c) has a multiple root}. (A. 5) 

The discriminant of P{a,c) is a polynomial in the parameters a — (ao,...,a„) z.e. Dp{a) — D{P) e 
IR[ao, . . . , a„], hence we can write 

V = {a^ (ao, . . . , a„) e R"+i : Dp{a) - 0}. (A.6) 

Such a set is, by definition, a proper Zariski-closed subset of R"+^. 

As the complement of a proper Zariski-closed subset of M"+^ is open and dense also in the Euclidean 
topology, then Hypothesis [3] is generic. ■ 



B Proof of Lemma [T7 



First we shall prove by induction on k that for all 9 G 6fc,o,/3o has 



< A/(fc-f),J- ( 1 + ^ ) , (B.l) 



for all k> -f 1. 

For k = + 1 one has 

P^o'"^'^^-^ E ^(j,mo,...,m,„,m)Q.„,cr...C (B.2) 

siP+mol)o + ---+mio f)io=Sio+l 

SO that any tree 9 contributing to /3q'''°^^' has si + 1 nodes and j leaves, hence \L{9)\ — si + 1 + j. 
Notice that q < [)o < t)i < . . . < f)^^ , hence one has 

\Li9)\ = l + si+j <l + '-^^. (B.3) 

q 

Moreover for k — t)ig + 1 the r.h.s. in (|B.ip is equal to 2 -f Si^/q, so that the bound (jB.ip holds, 
because one has q > 1. Assume now that the bound (jB.ip holds for all k' < k and let us show that then 
it holds also for k. 

We call Mq — Mi)ig -I- 1 + Si^/q, so that the inductive hypothesis can be written as 

\L{9')\ < Mk{9') - Mo, (B.4) 
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for all 9' G 0fc',o,,3o: < k- 

By the inductive hypothesis, we have (cf. Figure [5]) 



\L{e)\ < 1 + si + so ~ s'oMo + MJ2H^^), (B.5) 



1=1 



for suitable 9i, . . . ,9si^ depending on 9. 

Let us set m -.^ k — t)i„ > 1. Hence, via the conditions ()6.5p we can write (IB.5|) as 



\L{9)\ < 1 + si + So - s'„Mo + M \^s,„ + m - siP - ^ soai)^j ■ (B.6) 

Hence we shall prove that 

1 + si + So - s'qAIo + M I Si, + m - sip - V so,^^)^ 1 < mM - 1 - (B.7) 

V t'o J 1 

or, in other words 

( siP + V so.;i)i] M + s'oMo > 5,,M + so + si + ^ + 2, (B.8) 

for all So, Sq, si > admitted by conditions (|6.5p . 

First of all for sf, = by the first condition in (|6.5p we have sip + so,o()o + ■ ■ • + so,io[)io — Si, + m. 
Moreover (si + so)q < sip + so^i„\)o + . . . + so.iobio = ^ia + '^i hence 

si + so<^^^, (B.9) 

q 

so that one obtain (jB.8|l if 

mM>2^ + 2 + -, (B.IO) 

q q 

hence one needs 

m 12- 

q 



m ( 2^ + 3) > 2^ + 2 + -, (B.ll) 



that is satisfied for all 771 > 1. 

For Sq = 1 the first conditions (|6.5p can be written as sip + so.ofjo + . . . + (so,io + ^)^ia — ^io + s° 



that 

5i+5o<Ii2 (B.12) 

Hence we obtain (IB. 81) if 



nM > + " '^'^ + 1. (B.13) 

q 

and again (|B.13|1 is satisfied because n = k — ki > 1 . 

Finally for Sg > 2 the first condition in ()6.5p can be written sip + so,o[)o + • • ■ + (so,;,, + So)hio ^ ^io^ 
so that si + So < Si + Sq + Sq < and we obtain (jB.8p by requiring 

s^,M + s'J^ + l)> s,„M + 2^ + 2, (B.14) 

V q / q 
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that is satisfied as we are assuming Sq > 2. 

This exhausts the discussion over all the choices of sq, s'q, si. 
Let us show now that 

\L{e)\<Mk-l, (B.15) 

forall0eefe,,j,/ = /3,B, fc>p. 

Again recall that a tree 9 G G>k,uj contributes to /i,*^' with f = /3, B, so that the bound (jB.15p is 
trivially satisfied for k = p because one has \L{6)\ = 1. 

Let us suppose now that the bound holds for all p < fc' < fc; again we shall prove that then it holds 
also for k. 

Recall that a tree contributing to /i,'^' is of the form depicted in Figure [51 where so,a is the number 
of the lines exiting a leaf with leaf label a and entering Oq, so = + . . . + so,io' ■^i is the number of the 
lines exiting a node and entering Oq, and Sq, s'l are the graph elements entering Op with component label 
Po and / respectively. Hence, by the inductive hypothesis and by the bound (|B.1|) . we have 

\L{e)\ < 1 + So + si - SqMo -s[+M J2 ^(^0' (B-16) 

for suitable 0i, . . . ,9gi^ depending on 0. 

Let us supposte first bo„ = 1; thus, via the first condition in (|6.2p . we have to prove the bound 

1 + So + si + M(fc - P - so,of)o - ... - so,,of)^o - sip) - s^Mq - s[ < Mk - 1, (B.17) 

or, in other words, 

io 

soAMl)^ - 1) + si(Mp - 1) + Mp + sfjAfo + s'l > 2, (B.18) 

i=0 

and this is obviously satisfied as Mt)i, Mp > 3. 
Finally if 6oo = we have 

^fc(0,) = fc-sip, so,o + ... + so,»o+So = 0. si+s'i>2, (B.19) 

i=l 

SO that, by the second condition in (|6.2p . we have to prove the bound 

1 + si + M{k - Sip) - s'l < Mk ~ 1, (B.20) 
or, in other words si(A/p — 1) + s'^ > 2, and again this is obviously satisfied as si + s'^ > 2 and Mp > 1. 
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